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ABSTRACT 

Gravitational lens time delays depend on the Hubble constant, the observed 
image positions, and the surface mass density of the lens in the annulus between 
the images. Simple time delay lenses like PG1115+080, SBS1520+530, B1600+434, 
PKS1830-211 and HE2149-2745 have Hq = A(l - («)) + B{n){r} - 1) where the two 
coefficients A ~ 90 km/s Mpc and B ~ 10 km/s Mpc depend on the measured delays 
and the observed image positions, (k) is the mean surface density in the annulus 
between the images, and there is a small correction from the logarithmic slope ry ~ 2 
of the surface density profile, k oc R 1 " 11 , in the annulus. These 5 systems are very 
homogeneous, since for fixed Hq = WOh km/s Mpc they must have the same surface 
density, (n) = 1.11 — 1.22h ± 0.04, with an upper bound of a K < 0.07 on any dispersion 
in (k) beyond those due to the measurement errors. If the lenses have their expected 
dark halos, (k) ~ 0.5 and Hq ~ 51 ± 5 km/s Mpc, while if they have constant 
mass-to-light ratios, (k) — 0.1-0.2 and Hq ~ 73 ± 8 km/s Mpc. More complicated 
lenses with multiple components or strong perturbations from nearby clusters, like 
RXJ0911+0551 and Q0957+561, are easily recognized because they have significantly 
different coefficients. 



Subject headings: cosmology: gravitational lensing; cosmology: Hubble constant; dark 
matter 



1. Introduction 

Gravitational lens time delay measurements can determine the Hubble constant Hq given a 
model for the gravitational potential of the lens galaxy (Refsdal |1964| ) . With 9 well-measured time 



delays (see Schechter 2000 ), gravitational lenses have become a serious alternative to the local 



distance scale. Kochanek ( 2002a| ) modeled the 5 simple, well-observed time delay lenses to find 



that Hq = 48^4 km/s Mpc if the lenses have isothermal mass distributions (flat rotation curves) 
and that Hq = 71 ± 6 km/s Mpc if they have constant mass-to-light (M/L) ratios. Intermediate 
values can be found by adjusting the dark matter halo between these two extreme models, but 
the results agree with the local estimate of Hq = 72 ± 8 km/s Mpc by the HST Key Project 
(Freedman et al. 200 1|) only if the lenses contain little or no dark matter. Models of other lenses 
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(e.g. Munoz, Kochanek & Keeton 2001| ), local stellar dynamical measurements (e.g. Rix et 
al. [1997 , Romanowsky & Kochanek |1999| , Gerhard et al. 2001, Treu & Koopmans |2002j ), weak 
lensing (e.g. Guzik & Seljak |2002| ) and X-ray (e.g. Fabbiano |j~98g| , Lowenstein & White [T999J) 
measurements all require mass distributions close to the dark matter limit. 

The dependence of the estimate of Hq from a gravitational lens time delay on the mass 
distribution of the lens is well known, both from models of particular time delay lenses (e.g. Keeton 
& Kochanek 1997 , Impey et al. 1998 , Koopmans &; Fassnacht 1999] , Lehar et al. 200C| , Keeton et 
al. 200C , Williams & Saha 200C, Winn et al. |2002 ) and from general analytic principles (Falco 



et al. |1985| , Gorenstein et al. |1988| , Refsdal & Surdej [1994 Witt, Mao & Schechter |1995| , Witt, 
Mao & Keeton |2000| , Schechter |2000| , Rusin |200q , Saha |2000| , Zhao & Pronk |2001| , Wucknitz |2002 
Oguri et al. 2002; ). The most important trend is that the predicted time delays, or the inferred 
Hq, increase as the mass distribution becomes more centrally concentrated. In particular, Falco 



et al. (1985) showed that adding a constant surface density sheet to the mass distribution has no 
observable effect other than to rescale the time delay (the mass sheet degeneracy), and Witt et 
al. Q2000| ) showed that the Hubble constants estimated from lens potentials of the form <f) oc R 5 ~ v 
(p oc r _?? ) roughly scale as Hq oc (r/ — 1)/At, almost independent of the angular structure of the 
potential. 

Our objective in this paper is to clarify the physical properties of the lenses which control the 
expected time delays after making full use of the astrometric constraints available for a typical 
lens. We will show that time delays are controlled by a local property of the lens, the average 
surface density in the annulus between the lensed images. Reducing the model-dependence of 
the time delay estimates to such a simple physical property of the lens leads to simple, accurate 
scaling laws for Hq, provides model-independent tests for the homogeneity of the time delay lenses, 
demonstrates that standard parametric models have the necessary degrees of freedom to study 
degeneracies in estimates of Hq from time delay measurements, and makes it easy to theoretically 
predict time delays for standard halo models. In §2 we develop analytic results for lenses in 
arbitrary spherical density distributions. We first consider circular lens potentials (§2.1) and then 
examine the effects of adding an external shear (§2.2), the quadrupole of a singular isothermal 
ellipsoid (§2.3) or a general quadrupole (§2.4). These analytic results show that the time delays are 
largely determined by the mean surface density {k) = (£) /S c in the annulus between the images. 
In §3 we show that the Hq estimates for most lenses can be reduced to a linear function of {k), 
and that these simple scaling solutions reproduce the full numerical models by Kochanek ( |2002a| ) . 
In §4 we discuss the consequences of these results for the homogeneity of the lens population. 



2. Analytic Models 



We start from the dimensionless time delay (see Schneider, Ehlers Sz Falco 1992| ), 



1 2 1 2 

t = -(x-u) +-{y-v) -<j)(x,y) 



(1) 
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where x and y are the angular coordinates on the lens plane, u and v are angular coordinates on 
the source plane, and 4>(x, y) is the lens potential. The center of the lens galaxy is used as the 
coordinate origin. The true delay is 

1 + ZiDolDoS f. rr-l fr) s 

* = tfo r H D LS T = D °" H ° T (2) 

where the Dij are angular diameter distances between the observer, lens and source, z\ is the 
lens redshift, and H§ is the Hubble constant. [] The factor containing the distances and the 
Hubble radius, rn = c/Hq, depends on the cosmological model, the source redshift z s and the 
lens redshift Z[, but not the Hubble constant. We assume a Oo = 0.3 flat universe, but the 
cosmological dependence of the results is very weak. The lens potential satisfies the Poisson 
equation V 2 (f) = 2k where k = S/S c is the surface density of the lens in units of the critical 
surface density S c = c 2 Dos I^kGDolDls- The first two terms of the time delay (Eqn. |l]) form 
the geometric delay from the bending of the rays, and the remaining term is the Shapiro delay 
from passing through the gravitational potential of the lens. We observe images at the solutions 
of Vr = 0. 

We derive the time delay differences for a pair of images in a series of analytic models 



based on multipole expansions for the lens potential (see Kochanek |1991| and Trotter, Winn & 



Hewitt 2000 for other applications of multipole methods to gravitational lenses). Our objective 
is to achieve expressions for the delays which are accurate to ^ 5% when estimating the Hubble 
constant. Two images at radius R\ and R2 have an average radius of (R) = (R\ + R-z)/2 and 
bracket an annulus of width AR = R,\ — R%- For the monopole of the lens, our accuracy goal 
will require expansions of order (Ai?/(i?)) 3 . The ellipticity of the lens potential generally satisfies 
e < AR/(R) for two-image lenses and e ~ AR/{R) for four-image lenses. Assuming a standard 
multipole sequence for a roughly ellipsoidal mass distribution, the quadrupole (exp(2i#)) is of 
order e and can be treated to one lower order of expansion than the monopole, and higher order 
poles (exp(2mi#)) can be neglected since they are of order e m . We start in §2.1 by estimating the 
time delays for an arbitrary spherical lens. Next we add a quadrupole potential to the lens: in §2.2 
we consider an external (tidal) shear, in §2.3 we consider the quadrupole produced by a singular 
isothermal ellipsoid, and in §2.4 we consider a quadrupole with an arbitrary ratio between the 
internal and external contributions to the quadrupole near the images. 



We note in passing that the distance factor scaling of the time delays is simpler if comoving distances are used 
rather than angular diameter distances. The extra factor of 1 + zi is absorbed into the distances. Comoving distances, 
not angular diameter distances, are the natural cosmological distances for gravitational lensing calculations as they 
almost always lead to simpler analytic expressions. 
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2.1. Spherical Models 



We start with a lens having two images on opposite sides of a spherical lens, 4>(x,y) = (po(R), 
defined by a monopole potential 4>o(R). The outer image, at radius R\, is at a minimum of the 
time delay surface and the inner image, at radius R 2 , is at a saddle point, and the two images and 
the lens lie on a line. Starting from Eqn. ([[]), the geometric delay between the two images can be 
written 

Ar 12j9eom = -i (Rl - (l - 2<«>) (3) 

where 

is the average surface density in the annulus bounded by the images in units of the critical surface 
density. Only the derivative of the potential is needed to determine the geometric delay, and the 
mean surface density enters because the radial derivatives (f)' of the monopole at the positions of 
the two images are related by 



(5) 



independent of the density distribution in the annulus. For an isothermal lens (k) = 1/2 and the 
geometric contribution to the delay is identically zero for all images. The Shapiro delay difference 
is more complicated, 



AT 12 ,Shap = ~<«) (Rl ~ Rl) + RlR2 (l 



1 -^2 



r 2 



Ri 



udu (^k(i 



( K )jln— , (6) 



because the Shapiro delay depends on the surface density distribution in the annulus as well as 
the mean surface density. The first term, which depends only on (k), is the most important, as 
the contribution from the integral is smaller by (AR/ (R)) 2 . Combining the terms, the time delay 
between the images, 



An 2 = n -T 2 = (l-<«>) -\ (R 



Rl) +R 1 R 2 ln^ 
7 Ri. 



R 2 , 
2 / udu (k(u) 
IRi v 



depends only on the image positions and the surface density of the lens in the annulus between 
them (Gorenstein et al. 1988). The first term is the same as the time delay of a point mass lens 



combined with a uniform sheet of density (k). 

We can capture the essential elements of Eqn. (|7]) without the complexity by assuming that 
the surface density in the annulus between the images is a power law k = Ki(i?/i?i) 1_r? and then 
expanding the time delay as a series in the ratio of the thickness of the annulus AR = R 2 — Ri to 
its average radius (R) = (Ri + R 2 )/2 to find that 



Ar- 



12 



Ri 



r\ 



!-<«> 



t](k) [AR\ 

\JR)J 



12 



+ 



(8) 
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The first term is the exact delay for a singular isothermal sphere (SIS, rj = 2, (k) = 1/2), and the 
series converges rapidly even for large values of AR/ (R). For example, the first term approximates 
the time delay of a point mass lens with only a 9% error even as AR/(R) — > 1. The second order 
correction adjusts for the small changes in the critical radius (inside which the average surface 
density is unity) with the shape of the density profile. It comes only from the Shapiro delay 
(Eqn. 0) and can be decomposed (1 — t}(k) = (1 — (k)) + {k)(i] — 1)) into a contribution depending 
only on the average surface density and a contribution which depends on its distribution. The 
dependence on the logarithmic slope r\ arises only from the integral term of the Shapiro delay. 
The scaling observed by Witt et al. ( 2000 ) for potentials of the form <fi tx R 3 ~ v arises because 
1 — (k) = (rj — l)/2 near the critical line of the lens. 



2.2. A Spherical Lens in an External Shear 

Real lenses are not spherical, so we next consider a range of models with a quadrupole as well 
as a monopole. The monopole depends only on the total mass inside the Einstein ring and the 
surface density in the annulus. The quadrupole depends on the total internal quadrupole (due to 
the quadrupole moment of the mass inside the Einstein ring), the total external quadrupole (due 
to the quadrupole moment of the mass outside the ring) and the quadrupole component from the 
surface density inside the annulus. Because the ellipticity of the potential is small AR/(R)) 
we can treat the quadrupole to a lower order of expansion than the monopole and ignore the 
quadrupole component of the surface density in the annulus. For our analytic results, but not for 
the numerical results in §3, we assume that the internal and external quadrupoles are aligned. 
Thus, we need three numbers, two amplitudes and an orientation, to define the quadrupole. 
Higher order angular structures (octopole etc.) will be smaller than the quadrupole by one or 
more powers of the ellipticity of the potential, and can be neglected. 

We start with two familiar examples of quadrupoles. We first consider an external (tidal) 
shear, where the quadrupole is generated entirely by material outside the region with the images. 
In §2.3 we discuss the quadrupole produced by a singular isothermal ellipsoid. Then in §2.4 we 
discuss the general aligned quadrupole. Ideally, we would present simple results to the same 
order ((AR/ (R)) 3 ) as we did for the monopole. However, our objective is to illustrate the effects 
of quadrupole structures on time delays, so we will not reproduce terms which are grotesquely 
complicated simply to reach an arbitrary order of expansion. 

The quadrupole potential associated with an external shear is <fe = (lext/2)R 2 cos 2(6 — 6 ex t), 
where j ex t is the amplitude of the external shear and 6 ex t is its orientation. As in §2.1, the 
monopole is defined by the mass inside the Einstein ring, the average surface density (n) in the 
annulus and the logarithmic derivative of the surface density rj. We model two images located 
at (Ri cos 0i,Ri sin#i) and (i?2 cos #2,-^2 sin 62) from the lens center. Using the positions of two 
images we solve for the mass inside the Einstein ring and the amplitude of the external shear. 
This results in expressions which depend on the surface density in the annulus ((k), rj) and the 
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orientation of the shear (9 ex t). The shear axis seems the best choice for an independent variable in 
our analytic models because its value is little affected by changes in the monopole or quadrupole 
structure of the lens (see Kochanek |1991| ). If we again expand the delay as a series in ARJ (R), 



At 



12 



r »+w ri+ (w) n 



then the leading term, 
T = 



1 



K) sin 



A9 



12 



(l - (k)) as A9 12 -» it, 



(9) 



(10) 



changes from the result for a spherical system only in adding a dependence on the angle between 
the images A9\ 2 = 9\ — 02- The limit as A9± 2 — > vr shows that we recover the result for a circular 
lens and images collinear with the lens (Eqn. |8|). The first order correction, 



!(< 



1 - (k)) sin (A0ia) cot(0i + 6 2 - 29 



ext J 



as AO 



depends on the shear axis 9 ex t, but the second order correction, 



To 



1 

12 



1 



(4 + 3cosA0i 2 ) + (k) (r]-l] 



1 - v{k) 
12 



as 



12 



A0 



7T, 



12 



7T, 



(11) 



(12) 



depends only on the surface density in the annulus. Note that the higher order terms return to 
the spherical case for images collinear with the lens. The amplitude of the shear is 



lext 



[1 - <«» (Rl - Rl) 



sin A9\2 



{Rl + R\- 2R X R 2 cos A0 i2 ) sin {9 1 + 9 2 - 29. 



as A9 



12 



7T, 



(13) 



ext ) 



so the singularity in T\ when 9\ + 9 2 — 29 ex t — ► is due to a singularity in the shear that will not 
arise in practice since 9 ex t must be set using the parameters of a reasonable macro model. Also 
note that the amplitude of the shear scales as 1 — (k), so that lenses with lower surface densities 
require larger shears. 

The dependence of the delay on the structure of the surface density in the annulus is identical 
to that for the spherical lens, with the integral in the Shapiro delay (Eqn. ||) adding a term 
(ji — 1)(k) /12 to the T 2 coefficient. All other terms in the time delay depend only on the derivatives 
of the monopole and are functions only of (k) even for an arbitrary monopole surface density. The 
term vanishes for rj = 1 when k = K\(R/ Ri) 1 ^ 71 becomes a constant surface density. 



2.3. A Spherical Lens in an SIE Quadrupole 

Unlike our result for an external shear, where the time delay depends on the angle between 
the images even at the lowest order of the expansion, Witt et al. ( p000| ) found that the time delays 
for potentials of the form oc RJ-(9) (which includes the standard singular isotherma ellipsoid or 
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SIE potential) depend only on the image radii. These models have a quadrupole potential of the 
form 02 = cR(R) cos2(# — 6 e ). The leading term in the time delay, 

T = 1 - <«> (14) 

is independent of the angle, just as for the Witt et al. ( 200CQ potential. The second term is a 
correction for the effects of the ellipticity, 

Ti = -- (l - («}) [5 cos (A0i 2 /2) + cos (3A0i 2 /2)] cot (0i + # 2 - 20 eiri ) esc 3 (A0i 2 /2) (15) 

The next term is too complicated to present in full, but it has the structure 

T 2 = (l - («)) g(9i,e 2 ,e e ) + A fa - (16) 

The first factor depends on a very complicated angular function, and the second factor is the 
correction to the Shapiro delay from the logarithmic slope rj of the monopole surface density. To 
lowest order, the ellipticity of the potential is 

_ 1-( K )AR^ AQ , 



(R) 



sin A0 i2 csc 4 (A#i 2 /2) csc(0i + 2 - 26 e ). (17) 



We again find that the ellipticity scales as 1 — (k) and that the singularities in T\ correspond to 
singularities in e. 



2.4. A Spherical Lens in a General Quadrupole 

The quadrupoles of an external shear and an SIE differ in the balance between internal and 
external quadrupoles. Since they also produce different angular dependences for the leading term 
in the time delay, it suggests that the angular dependence of the time delay is determined by 
the relative strength of the two shear components. We can test this by determining the leading 
term in the time delay for a quadrupole which is an arbitrary sum of an internal quadrupole of 
amplitude and external quadrupole of amplitude "fext, 

$2 = \ [lextR 2 + lint^\ COS 1(0 - 7 ). (18) 

If we use the image positions to determine the total shear, T = ji n t + 7ezt, as a function of the 
shear angle # 7 and the ratio fi n t = "fint/T, we h n d that the leading term, 

n sin 2 (Afl 12 /2) 
{K) ) 1 - 4/ mi cos* (A0 12 /2) ' Uyj 



has the standard scaling with 1 — (k) and an angular dependence determined by the balance 
between the internal and external quadrupoles. For an external shear (fi n t = 0) the denominator is 
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unity, and for the quadupole of an SIE (fi n t = 1/4 for ki oc R 1 ) the numerator and denominator 
cancel to make To independent of A9\2 just as in the Witt et al. p000|) results. For f int > 1/4, 
where the internal shear begins to dominate, there are angles where the expansion fails. The 
higher order terms are too complex to be worth presenting for general configurations. However, 
for two-image lenses where the images lie on opposite sides of the lens, 62 = 0\ + it + 56, with a 
small angular offset defined by £ = 58(R)/AR <S 1, the delay depends little on the structure of the 
quadrupole. The expansion terms, 

T = 1-(k), T 1 = and T 2 = ~(1 - v ){k) + ^ (l - («)) [2cos20 7 +£(4/ int - 1)] (20) 

are nearly identical to those for a spherical lens (Eqn. |||) because the quadrupole structure enters 
only at second order. 



3. A Simple Semi- Analytic Model 

These analytic results suggest a simple semi-analytic model for understanding the time 
delays observed in gravitational lenses. It cannot be a completely analytic model because we 
require a numerical solution using all the lensed images to determine the quadrupole structure 
of the four-image lenses. We continue to model the monopole by the mass inside the average 
image radius and a power law, k(R) oc ^ for the surface density between the inner and outer 
images parameterized by the mean surface density (k) and the logarithmic slope r/. The resulting 
expression for the monopole deflection, 

(R) 2 r R 

= bo V + p / <R)RdR, (21) 
R R J(r) 

where 60 ~ (R) determines the mass inside (R) is very similar to an early model of Q0957+561 by 
Borgeest & Refsdal ( |1984| ) which combined a point mass and an isothermal sphere (77 = 2). Since 
many of the lenses have AR/(R) ~ 0.5 we expanded the monopole to first order in the deviation 
of the density exponent from isothermal, 5rj = rj — 2, rather than AR/{R). This leads to results 
that are exact for rj = 2 independent of the width of the annulus, and are reasonably accurate for 
1^??^3. We combined the monopole with a general quadrupole 

MR, 0) = \l ex tR 2 cos 2{9 - 9 ext ) + \lint^r cos 2(6 - 6 int ) (22) 

where the internal and external components have amplitudes and orientations of (jtntfiint) an d 
{lextfiext) respectively. For each lens we determined the values of bo and the shear components 
which minimized the differences in the projected source positions, leaving expressions depending 
only on the monopole structure in the annulus ((k) and rj). For the current time delay lenses, we 
cannot determine (k) or r\ from the available constraints. By doing the fits on the source plane, 
the solutions can be found rapidly using mathematica programs. We did not fit the image flux 
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ratios. For the two-image lenses, where there are too few constraints to determine all 5 variables, 
we stabilized the solutions by finding the minimum shear solutions.^] We experimented with the 
structure of the shear terms (and the shear prior for the two-image lenses) to ensure that the 
results reported for Hq are not affected by the assumed quadrupole structure. For the two-image 
lenses the results had no significant dependence on the quadrupole structure (as expected from 
§2.4), and for the four-image lenses the quadrupole structure was well-determined by the image 
positions. Errors were estimated by Monte Carlo simulations and for Hq they are dominated by 
the errors in the time delay measurements. 

Once we have solved for the model parameters required to match the image positions, we 
derive the time delays between images as a function of the remaining parameters (re) and r\. Based 
on the analytic and numerical results, we can express the results as simple scaling relations for the 
Hubble constant of the form 

H ~A(l-{K))+B{ K )(r)-l) + C (23) 

where the three coefficients (A, B and C) have the units of km/s Mpc and scale inversely with the 
measured time delays. We have isolated the terms more physically here than in the expansions of 
§2, isolating the dependence of the logarithmic slope r\ in the B coefficient, but the relationship 
between the coefficients and the expansion can be worked out for each model. For example, the To 
term defined in Eqn. @) contributes — (R\ — B%)TQD e ff / At\2 to the A coefficient, where At\2 is 
the measured delay and D e ff is the effective distance defined in Eqn. (||). In the analytic models 
C = 0. For simple lenses the Hubble constant Hq ~ A for very centrally concentrated lenses 
((k) — > 0), and Hq ~ (A + B)/2 for an isothermal model ((re) = 1/2, r/ = 2). We order the images 
so that the time delay is positive and A > 0. The effects of an additional mass sheet of density n ex t 
can be easily added to the models. The adjustments depend on whether we regard (re) as a fixed 
number or a fit parameter whose value is also changed by the addition of the extra convergence. 
If (re) is a fixed number, then the estimate of the Hubble constant changes to 

Hq ~ (1 - K ext ) [A (l - (re)) + C] + B(k)( V - 1). (24) 

In most cases, however, we start with a lens whose mass and surface density are scaled to 
fully explain the image splitting. With the addition of an external convergence, the annular 
surface density (k)o of this reference model is reduced to (re) = (1 — K ext )(K)o because the extra 
convergence reduces the mass of the primary lens. Thus, if we specify (re)o, the estimate of the 
Hubble constant with the addition of an external convergence is 

H ~ (1 - K ext ) [A (l - (re)o(l - K ext )) + 5(re) (r? - 1) + c] (25) 

and it has the familiar Hq oc (1 — K ex t) scaling of the mass-sheet degeneracy (Falco et al. 1985| ). 



2 For two-image lenses we added on weak prior on the shear amplitudes to the \ 2 for the mapping 
of the images back to a common source. The prior was defined by T ■ Tj^ 2 where T = 
{jint cos 28int, 7int sin 2#i„t , j e xt cos 28 ex t, Jext sin 20 e xt) is the vector of Cartesian shear components and we set 
70 = 0.05. 
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Table 1. Hubble Constant Models 



Lens 


Images 


Delay 




Re/(R) 


AR 

m> 


A 


A 


B 


C 






(days) 




[Ho] 


[Ho] 


[Ho] 


[H ] 


RXJ0911+0551 


AD 


146 ± 


4 


0.55 


0.58 


272 ± 14 


82 ± 5 


69 ±4 


-9 ± 1 




BD 








0.54 


280 ± 14 


81 ± 5 


69 ±5 


-9 ± 1 




CD 








0.60 


262 ± 14 


82 ± 6 


68 ±5 


-9 ± 1 


Q0957+561 




417 ± 


3 


0.63 


0.67 


249 ± 2 


221 ± 2 


51 ± 1 


-7±0 


PG1115+080 


A]B 


12 ± 


2 


0.41 


0.10 


81 ± 17 


70 ± 14 


0±0 


1 ±0 




A X C 


13 ± 


2 




0.10 


126 ± 22 


112 ±20 


9±2 


0±0 




BC 


25 ± 


2 




0.19 


97 ± 8 


93 ± 8 


5±0 


0±0 


SBS1520+530 




129 ± 


3 


0.51 


0.52 


98 ± 5 


92 ± 5 


10 ±1 


-1 ±0 


B1600+434 




51 ± 


2 


0.94 


0.64 


119 ±14 


104 ± 10 


21 ±9 


-3 ± 2 


PKS1830-211 




26 ± 


4 




0.19 


85 ± 15 


88 ± 16 


1±0 


0±0 


HE2149-2745 




103 ± 12 


0.55 


0.59 


95 ± 12 


84 ± 11 


14 ±2 


-1 ±0 



Note. - For each lens and image pair we give the measured time delay in days, followed 
by the coefficients which determine the Hubble constant given the surface density model. We 
first give the single coefficient Aq needed for the lowest order expansion Hq = Aq(1 — (k)) for 
an arbitrary monopole in an external shear, and then the three coefficients A, B and C used 
in our standard expansion (Eqn. [23] ) for an arbitrary monopole and quadropole. The formal 
errors include the uncertainties in both the astrometry and time delays but are dominated by 
the time delays. This means that the uncertainties are almost perfectly correlated! We broaden 
the time delay errors to a minimum error of 5% to account for systematic uncertainties such 
as convergence fluctuations from large scale structure. The ratio R e /(R) is the ratio of the 
(monopole) half-light radius to the average image radius, and the ratio AR/2(R) is an estimate 
of the thickness of the annulus. 



Table 2. Comparisons to Full Numerical Models 



Lens Image Full Model Hq {K)o/n ex t Lowest Order Hq Standard Hq 







DM 


M/L 


DM 


M/L 


DM 




M/L 




DM 


M/L 




RXJ0911+0551 


AD 


49 


±4 


67 ±4 


0.50/0.24 


0.11/0.31 


129 ± 


7 


174 ± 


9 


58 ± 


3 








BD 








0.50/0.24 


0.10/0.31 


132 ± 


7 


180 ± 


9 


58 ± 


3 








CD 








0.50/0.24 


0.12/0.31 


123 ± 


7 


166 ± 


9 


58 ± 


3 






PG1115+080 


AiB 


46 


±4 


69 ±5 


0.50/0.09 


0.18/0.16 


40 ± 


8 


58 ± 12 


35 ± 


7 


51 ± 10 




A X C 








0.50/0.09 


0.11/0.16 


62 ± 11 


96 ±17 


60 ± 


10 


87 ± 15 




BC 








0.50/0.09 


0.15/0.16 


48 ± 


4 


71 ± 


6 


48 ± 


4 


69 ± 


6 


SBS1520+530 




53 


±6 


73 ±5 


0.50/0.00 


0.20/0.00 


49 ± 


3 


79 ± 


4 


51 ± 


3 


77 ± 


4 


B1600+434 




59 


±6 


71 ±9 


0.50/0.05 


0.43/0.18 


60 ± 


7 


63 ± 


7 


59 ± 


7 


68 ± 


9 


HE2 149-2745 




47 


±8 


66 ±8 


0.50/0.00 


0.22/0.00 


48 ± 


6 


74 ± 


9 


48 ± 


6 


70 ± 


9 



Note. - - Hubble constant estimates for dark matter (DM) and constant mass-to-light ratio (M/L) models 
of the density distribution. All estimates of the Hubble constant are in units of km/s Mpc. The estimates from 
full numerical models are the dark matter and constant M/L models from Kochanek ( [2002a ). The lowest order 
Hq estimates use the lowest order expansion Hq = Aq(1 — («)), while the standard Hq estimates use Eqn. (j23|). 
The (re)o/ftext columns give the mean surface density (k)q of an unperturbed spherical model and the surface 
density due to external perturbers K ex t in the full numerical models. For the isothermal models we used rj = 2 
while for the constant M/L models we used rj = 3. 
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We computed these simple models for 7 of the 9 time delay lenses, excluding only B0218+357 
(where the lens position is poorly known, Lehar et al. |2000| ) and B1608+656 (where the lens 
consists of two interacting galaxies, Koopmans & Fassnacht |1999| ). We consider RXJ0911+0551 
(Hjorth et al. [200% Q0957+561 (Schild & Thomson |l~997| , Kundic et al. [1997| , Bernstein et 
al. |l~997l , Keeton et al. p000| ), PG1115+080 (Schechter et al. $99% , Barkana [1997| , Impey et 
al. |i~998D , SBS1520+530 (Burud |2002b| , Faure et al. [20021) , B1600+434 (Burud et al. |2500| , 
Koopmans et al. |2000|) , PKS1830-211 (Lovell et al. |1998| , Winn et al. |2002| , Courbin et al. ( |2002D ) 



and HE2149-2745 (Burud et al. 2002 ). The system parameters are generally derived from the 



CfA/Arizona Space Telescope Lens Survey (CASTLES, Falco et al. 2001 ) images of the systems. 



Of these 7 systems, we expect our simple models to have difficulty with Q0957+561 and 
RXJ0911+0551. Q0957+561 is near the center of a cluster (Keeton et al. p000| , Chartas et al. |20"02"1) 
where an external shear will poorly describe the effect of the cluster on the image positions. 
RXJ0911+0551 is on the outskirts of a cluster (see Morgan et al. 2001| ) and the primary lens 
also has a small satellite. The external shear approximation is a poor one for the perturbations 
from the RXJ0911+0551 cluster. It is at best a marginal approximation for {k) ~ 1/2 and it fails 
completely for (n) 



(recall that 7 oc (1 — (k)), see Eqns. 13 and [L7p. We included an SIS satellite 
in our model of RXJ0911+0551, so it will match the Kochanek ( |2002a| ) dark matter models, but 
it is difficult to treat the two potentials self-consistently within our approximation. We also note 
that our estimates for PKS1830-211 are invalid if we adopt the Courbin et al. ( [2002] ) multiple-lens 



interpretation of the system rather than the Winn et al. ( 2002 ) single- lens interpretation. 



We now check our approximations against full numerical models for RXJ0911+0551, 
PG1115+080, SBS1520+530, B1600+434 and HE2149-2745 by Kochanek ( [2002^ . We first 
checked the dependence of the results on the assumed quadrupole structure by fitting models with 
either a pure external shear, an SIE quadrupole or a pure internal shear. For the two-image lenses, 
the Hq estimates are essentially independent of the quadrupole structure, as expected from the 
analysis in §2.4. The largest fractional changes in Hq are approximately 3%. This is fortunate 
because the observational constraints on the two-image lenses cannot distinguish between the 
quadrupole structures. When we allow a general quadrupole, which matches realistic models 
where there are contributions to the angular structure from both tidal shears and the ellipticity of 
the lens galaxy (Keeton, Kochanek & Seljak 1997] ), all the doubles can be fit perfectly and have 
Hq estimates consistent with the three constrained shear models. 

For the two four-image lenses, RXJ0911+0551 and PG1115+080, the quadrupole structure is 
enormously important because the functional dependence of the delay on the angle between the 
images is so strong once the images do not lie on opposite sides of the lens. For example, for an 
SIS monopole ((k) = 1/2, r] = 2), the Hq estimates from the BC image pair in PG1115+080 are 
47 and 67 km/s Mpc for an external shear and an SIE quadrupole respectively. The ratio of the 
values matches that expected from Eqn. (^). Because they are four- image lenses, however, they 
also strongly distinguish between quadrupole structures, with both lenses being far more consistent 
with an external shear than either a SIE quadrupole or an internal shear. Time delay ratios 
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are also controlled by the quadrupole structure, and the measured delay ratios in PG1115+080 
are also more consistent with an external shear than either an SIE quadrupole or an internal 
shear. The delay ratio seemed to be model independent in our studies of PG1115+080 (Keeton 



& Kochanek 1997, Impey et al. 1998]) because the models focused on changes in the monopole 



while always optimizing the quadrupole structure. Since the best fit quadrupole was close to an 
external shear no matter what we used for the monopole structure, the time delay ratios showed 
almost no model dependence. For models with a general quadrupole, we find that fi n t ~ 0.02 and 
0.01 for RXJ0911+0551 and PG1115+080 respectively. While very close to a pure external shear 
(fint = 0), the general quadrupole significantly improves the fit to the image positions. 

We adopt the general quadrupole models as our standard. For the two-image lenses the 
choice has negligible effect on the estimates of Hq, and for the four- image systems it is required 
to provide an acceptable fit to the image positions. Table 1 presents the coefficients A, B, and C 
(Eqn. |23| ) as well as the coefficient Aq, with Hq = Aq(1 — («)), we would obtain using the lowest 
order term for an external shear model (To from Eqn. |l0|). The uncertainties in the coefficients 
are almost entirely due to the time delay estimates, although they do include the astrometric 
uncertainties in the image and lens positions. These become important only for B1600+434 where 
the lens position is relatively uncertain due to the dust lane in the galaxy (Koopmans, de Bruyn 
& Jackson 1998 , Mailer et al. 2000) ), The dominant role of the time delay uncertainties means that 



the uncertainties in the coefficients are almost perfectly correlated. With the exception of the two 
lenses in clusters, RXJ091 1+0551 and Q0957+561, the lowest order coefficient Aq matches the 
higher order estimate A to approximately 10%, the B coefficient is needed to achieve 10% accuracy, 
and the C coefficient is negligible compared to other sources of uncertainty. The amplitude of B 
is largest for the lenses with thicker annuli (larger AR/(R)). In fact, the coefficients clearly divide 
the sample into a set of five very similar, relatively isolated lenses (PG1115+080, SBS1520+530, 
B1600+434, PKS1830-211, and HE2149-2745) and the two lenses in clusters (RXJ0911+0551 and 
Q0957+561). 

The real question, however, is whether these simple semi-analytic representations can 
reproduce the results of full numerical models. Table 2 compares the semi-analytic estimates to 
the full numerical lens models from Kochanek ( |2002a ) for either a dark matter dominated, SIE 



lens model ({k)o = 1/2, r/ = 2) or a constant M/L model where {k)q ~ 0.1-0.2 for most lenses. We 
defined (k)q for the constant M/L models by the average density in the annulus for a spherical 
model using the intermediate axis effective radius (the geometric mean of the major and minor 
axes) for the scale length in a de Vaucouleurs profile. For a lens with surface density profile T,(R), 
the critical surface density profile is simply 

2/,f c RT,(R)dR 

where R c is the average critical radius of the lens. It depends only on the shape of the profile 
because for a spherical lens the mean density inside the critical radius is equal to the critical 
density so the image geometry determines the mass scale. We also include the effects of the 
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external convergences n ex t from the numerical models. We do not include the constant M/L 
semi-analytic estimates for RXJ0911+0551 where the treatment of the nearby cluster as an 
external shear becomes untenable. 

First, note that the estimates for each image pair in the four-image lenses RXJ0911+0551 
and PG1115+080 are mutually consistent. This holds for the A 2 B and A 2 C pairs of PG1115+080 
as well. Second, for the simple lenses, the lowest order and standard approximations are nearly 
identical for the dark matter (high (k)) models, but the lowest order approximation tends to give 
higher estimates for Hq for the constant M/L (low (n)) models. For the more complicated lens 
RXJ0911+0551, the differences are enormous. Finally, the approximate solutions are an excellent 
match to the full numerical models using ellipsoidal lenses in external shear fields. Using the AD 
image pair from RXJ0911+0551 and the BC image pair from PG1115+080, the average difference 
between the approximate and numerical dark matter models is AHq = 2 ± 4 km/s Mpc with 
almost all the differences due to the most complex system, RXJ0911+0551. The approximate 
model for RXJ0911+0551 completely fails in the constant M/L limit (see above), but if we 
compare the remaining four lenses the difference is AHq = 2 ± 1 km/s Mpc. For comparison, 
the average differences for the lowest order approximation are AHq = ± 2 km/s Mpc and 
AHq = 2 ± 6 km/s Mpc respectively, but we have to drop RXJ0911+0551 for the dark matter 
models as well as the constant M/L models. Some of these residuals may be due to the ambiguities 
in defining (n) for ellipsoidal lens galaxies. We conclude that our approximate solutions generally 
provide estimates for Hq to better than 10% individually and to better than 5% as an ensemble. 



Discussion 



Gravitational lens time delays are determined by the Hubble constant, the positions of the 
lensed images, and the surface density in the annulus bounded by the images. The average surface 
density (k) in the annulus is more important than its distribution. The relationship between time 
delays and the local surface density is exact for circular lenses (Gorenstein et al. 1988] ). While it 



is not exactly true for non-circular lenses, it is true in practice. In two-image lenses, where the 
images lie on opposite sides of the lens galaxy, the delays are insensitive to the angular structure 
of the lens. In four- image lenses the delays are very sensitive to the quadrupole structure of 
the potential, but the image positions tightly constrain the quadrupole and leave the surface 
density as the only important variable. Our local interpretation of time delays agrees with earlier 
understandings based on global models and degeneracies, but better isolates the relevant physical 
properties of lens models. It also allows us to find simple scaling solutions for the Hubble constant 
which accurately reproduce the results of full numerical models. These scaling solutions can be 
used to classify time delay lenses, since their structure differs for simple lenses and lenses strongly 
perturbed by clusters, and to quickly estimate the Hubble constant predicted for any assumed 
radial mass distribution. 



We can also use the scaling solutions to study the homogeneity of the lens population. 
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Because Ho must be a universal constant, we can use the requirement that the lens all produce the 
same Hubble constant to estimate the surface density differences between the lenses. If we assume 
a universal density slope 77, then the five simple lenses (PG1115+080, SBS1520+530, B1600+434, 
PKS1830-211 and HE2149-2745) are consistent with a common value for the average surface 
density of 

(k) ~ 1 - 1.07/i + 0.14(77 - 1)(1 -h)± 0.04, (27) 

where H = lOO/i km/s Mpc. Formally, only limited ranges for the Hubble constant are consistent 
with a universal surface density, with h < 0.7, h = 0.5 ± 0.6 and h = 0.9 ± 0.3 for 77 = 1 (constant 
surface density), 2 (isothermal) and 3 respectively. An alternative way to characterize the 
similarities is use one lens as a reference and then estimate the surface density differences between 
the reference lens and the other systems. For example, if we assume 77 = 2 and use PG1115+080 
as the reference system, we find that (k) = 1.05 — 1.13/i ± 0.04 for PG1115+080. The surface 
density differences are Aft = 0.07 - 0.09/t ± 0.06, 0.19 - 0.07/i ± 0.11, 0.00 - 0.00/t ± 0.19 and 
0.14 - 0.28/j ± 0.15 for SBS1520+530, B1600+434, PKS1830-211 and HE2149-2745 respectively. 
The images in PG1115+080, SBS1520+530 and HE2149-2745 all lie at approximately the same 
radius relative to the lens galaxy (R e /(R) ~ 0.5) while the images in B1600+434 are close to the 
half-light radius (R e /{R) ~ 1). Reasonable slope differences between the lenses can change the 
surface density differences by 0.00-0.27(1 — h)Ar] depending on the lens. At least for constant 77, 
the upper bound on the existence of any scatter in the surface density beyond that implied by 
the measurement errors is a K ^ 0.07, so these five simple lenses have very homogeneous intrinsic 
properties. 

Our analysis is not well suited to the two cluster lenses, RXJ091 1+0551 and Q0957+561, 
because our use of an external shear to represent perturbations is breaking down. If, however, 
we assume that the true surface density of the primary lens is the same as for the other lenses, 
we can estimate the external convergence n ex t required to make the time delay estimates for the 
cluster lenses agree with the other 5 systems. We find that K ea; t,09ii — 1H ~~ 1.437i ± 0.10 and 
^ext,0957 — 0.99 — 0.45/i ± 0.01. While these results should not be taken at face value because 
of the problems with our approximation for lenses which are too close to the central regions of 
clusters, the pattern of the solutions seems reasonable. The RXJ0911+0551 lens lies 200/1" 1 kpc 
from the cluster center, while the Q0957+561 lens lies only 30/i _1 kpc from the cluster center 
so we would expect larger surface densities in Q0957+561 than in RXJ0911+0551. Presumably 
similar analyses with more realistic models for the cluster would improve matters. 

Our semi-analytic models for the four simple lenses with good photometry, PG1 115+080, 
SBS1520+530, B1600+434 and HE2149-2745, almost exactly reproduce the H estimates based 
on full numerical models by Kochanek ( 2002a| ). We find Hq ~ 51 ± 5 km/s Mpc if k = 1/2 



and 77 = 2 as expected in dark matter dominated lenses and Hq ~ 73 ± 8 km/s Mpc if they 
have constant mass-to-light ratios. The dependence of the time delays on the radial and angular 
structure of the lens shows that standard parametric models, which can adjust both the radial 
mass distribution and the quadrupole structure of the lens, encompass the physical properties 
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needed to study the dependence of Hubble constant estimates on the mass distribution of the lens. 
Most degrees of freedom in non-parametric models (e.g. Williams & Saha 2000| ) are not important 



for time delay estimates in simple lenses. Since the parametric models have the advantage of 
corresponding to physical models of galaxies, while most realizations of the non-parametric models 
do not (see Munoz et al. |2001| ), the use of non-parametric approaches is probably better suited 
to very complicated systems like B 1608+656 where the lens consists of two interacting galaxies 
(Koopmans & Fassnacht 1999| ) and it is unclear how to properly parameterize the system. It is 



also easy to make theoretical estimates for the time delays expected for standard cold dark matter 
(CDM) halo models, since we need only estimate the surface density (k) of the halo. We examine 
this problem in Kochanek ( 2002b|) . 
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